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(57) ABSTRACT 

A system and method of determining locations of lightning 
strikes has been described. The system includes multiple 
receivers located around an area of interest, such as a space 
center or airport. Each receiver monitors both sound and 
electric fields. The detection of an electric field pulse and a 
sound wave are used to calculate an area around each 
receiver in which the lighting is detected. A processor is 
coupled to the receivers to accurately determine the location 
of the lighting strike. The processor can manipulate the 
receiver data to compensate for environmental variables 
such as wind, temperature, and humidity. Further, each 
receiver processor can discriminate between distant and 
local lightning strikes. 

20 Claims, 1 Drawing Sheet 
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SYSTEM AND METHOD OF LOCATING 
LIGHTNING STRIKES 

This patent application is relate’! to U.S. Provisional 
Patent Application Ser. No. 60/182,404, entitled “Method 
and Apparatus for Accurate Location of Lightning Strikes , 
filed on Feb. 14, 2000. 

ORIGIN OF THE INVENTION 

The invention described herein w as made in the perfor- 
mance of work under a NASA contract and is subject to the 
provisions of Section 305 of the National Aeronautics and 
Space Act of 1958, as amended. Public Law 85-568 (72 Stat. 
435; 42 U.S.C. §2457). 

BACKGROUND OF THE INVENTION 

The present invention relates generally to identifying 
locations of lightning strikes. 

Electronic equipment is susceptible to damage caused by 
nearby lightning strikes. The accurate knowledge of a light- 
ning striking point is important to determine which equip- 
ment or system needs to be tested following a lightning 
strike. Existing lightning location systems can provide cov- 
erage of a wide area. For example, a lightning location 
system can provide coverage of an area having a 30 km 
radius. This system, however, has a 50% confidence region 
of about 500 meters. That is, the system has a 50% confi- 
dence that a lighting strike is within 500 meters of an 
identified location. As such, present lightning location sys- 
tems cannot be used to determine whether a lightning strike 
occurred inside or outside of a parameter of an area of 
concern. One such application of a lightning location system 
is a space shuttle launch pad for th< National Aeronautics 
and Space Administration (NASA). By accurately determin- 
ing lightning strike locations, electronic equipment located 
within the launch pad area can be test ed and/or reset to avoid 
erroneous operation. 

One method of determining the location of lightning 
strikes uses a set of video cameras that are pointed in 
different directions within the area ol concern. If a lightning 
strike occurs within the field of view of three or more 
cameras, the location of the strike can be determined. 
However, if the cameras are not pointed in the correct 
direction, or either an object or a heavy rain downpour 
obscures their field of view, it is difficult or impossible to 
accurately determine a striking point of the lightning. 
Further, this method has a relatively large uncertainty and 
does not facilitate an accurate location of the exact point of 
contact to the ground. 

For the reasons stated above, and for other reasons stated 
below which will become apparent to those skilled in the art 
upon reading and understanding the present specification, 
there is a need in the art for the system and method to 
accurately locate locations of lightmng strikes. 

SUMMARY OF THE INVENTION 

The above-mentioned problems with lightning strike 
location and other problems are addressed by the present 
invention and will be understood b> reading and studying 
the following specification. 

In one embodiment, a system to determine a location of 
lighting strikes comprises a processor, and a plurality of 
receivers coupled to the processor. Each of the receivers 
comprises an electric field sensor, an acoustic sensor, and a 
processor to provide a receiver output indicating a calculated 
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time differential between an electric field pulse and a sound 
wave (thunder). The processor determines the location of 
lighting strikes in response to the output from the plurality 
of receivers. 

s A method is provided for determining a location of 
lightning strikes. The method comprises locating a network 
of at least three electric field sensors and at least three sonic 
sensors in an area of interest, and collecting lightning strike 
information, including a difference of the time of arrival of 
o an electric field pulse and an associated sound w ave from a 
lightning strike. The method processes lightning strike infor- 
mation recorded by the at least three electric field sensors 
and at least three sonic sensors, including measuring a time 
difference between the arrival of the electric field pulse and 
5 the sound wave at each electric field and sonic sensor. The 
processor uses the time differentials to produce estimates of 
the range between the receiver and the lightning strike. The 
processed information is used to determine the location of 
the lightning strike. 

!() 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of one embodiment of a receiver 
of the present invention; 

>s FIG. 2 illustrates circles defined by a time difference 
between the arrival of an electric field signal and a sonic 
signal; and 

FIG. 3 is a block diagram of a system of the present 
invention, 
to 

DETAILED DESCRIPTION OF THE 
INVENTION 

In the following detailed description of the preferred 
embodiments, reference is made to the accompanying draw- 
ings that form a part hereof, and in which is shown by way 
of illustration specific preferred embodiments in which the 
inventions may be practiced. These embodiments are 
described in sufficient detail to enable those skilled in the art 
to practice the invention, and it is to be understood that other 
40 embodiments may be utilized and that logical, mechanical 
and electrical changes may be made w ithout departing from 
the spirit and scope of the present invention. The following 
detailed description is, therefore, not to be taken in a limiting 
sense, and the scope of the present invention is defined only 
by the claims. 

The present invention provides a system and method that 
can accurately locate a lightning strike within an area of 
interest. In one embodiment, the present invention can 
S() accurately locate a lightning strike within a few meters. 
Further, receivers used to detect lightning strikes can be 
located at distances of approximately one kilometer and 
greater apart. As explained below 7 , the present invention uses 
a combination of electric field and sonic sensors. 

55 The fast varying electric current associated with lightning 
discharges generates large electric field variations. The elec- 
tric field waveform propagates at the speed of light in a 
radial direction from the striking point of lightning. The 
sudden heating of the air caused by the large currents 
60 associated with the lightning discharge produces a sudden 
expansion of the air near the lightning channel. This results 
in a sound wave (thunder) that initially, for the first few 
meters, propagates at a supersonic speed and later propa- 
gates at a sonic speed. 

65 For an observer located remotely from a lightning strike 
location, the electric field waveform arrives earlier than the 
sonic sounds. This is because the electric field waveform 
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travels and a speed of approximately 300,000,000 m/s, while 
the sound wave travels and approximately 350 m/s. The 
observer can estimate the distance to striking point by 
measuring the time difference between the arrival of the 
electric field waveform and the arrival of the sound wave. 
This measurement defines a circle, with the observer at the 
center, on which the lightning strike might have occurred. 
The second observer and a different location using the same 
type of measurement also has a circle defined around them 
in which the lightning might have occurred. These two 
circles intersect at two points. With the addition of a third 
observer, a single striking point can be determined. The 
present invention provides receivers that can be located 
remotely from each other to accurately determine lightning 
strike locations. 

Referring to FIG. 1, a block diagram of one embodiment 
of a receiver 100 of the present invention is illustrated. The 
receiver includes an electric field antenna 110 , a microphone 
120 and a microcontroller 130 coupled to receive amplified 
signals from the antenna and microphone. In operation, 
multiple receivers are located around an area of interest. 
Each receiver is used to detect changes in an electric field 
and detection of sonic sounds. The combination of the 
electric field and sonic level are used to determine a location 
of the lightning strike. In operation, the electric field antenna 
110 detects changes in the electric field surrounding the 
antenna. The microphone 120 is used to detect sonic sounds 
surrounding the antenna. The output of each of these com- 
ponents can be amplified by amplifiers 140 and 150, if 
necessary, and the microcontroller 130 performs an analysis 
of the outputs to determine a radius around the receiver in 
which a lightning strike may have occurred. An example of 
this analysis is described in more detail below. By using 
three remotely located receivers, an accurate location of. 
lightning strikes can be determined. 

Referring to FIG. 2, the location of a lightning strike can 
be determined by the intersection of three circles 180, 190 
and 200 defined by the time difference between the arrival . 
of an electric field signal and a sonic signal. Each of the three 
circles are located around three remotely spaced receivers, 
210, 220 and 230, respectively. As explained above, each of 
the receivers determines a radius within which the lightning 
strike may have occurred. This radius, or distance d„, can be 
defined as the speed of sound (c) times the time difference 
(t„) between electric field signal in the sonic signal (d,= 45 

c*y. 

The present system relies on the accurate determination of 
the time elapsed between the reception of the electric field 
waveform and the reception of the sound waveform. The 
electric field waveform has a rise time in the order of a few 50 
microseconds, and its start can be determined with an 
accuracy of a fraction of a microsecond. At close distances, 
the sound waveform has a sharp wavefront. That is, its high 
frequency content (frequency>10 kHz) is a large. The sound 
waves propagating through air suffer large attenuation at 55 
high frequencies as compared to the attenuation of the low 
frequencies. The attenuation of the high frequencies is 
further enhanced by heavy rain. Thus, at close distances, the 
start of the sound wavefront can be easily measured since it 
has a fast rise time. At distances of greater than one or two 60 
kilometers, the wavefront rises slower, making it difficult to 
detect the exact time of the sound wave. 

In one embodiment of the present invention, the system 
includes a network of at least three receivers at different 
locations within a perimeter of interest to be monitored for 65 
lightning strikes. The microcontroller of each receiving 
station is used to measure the time difference between the 
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arrival of the electric field pulse and the arrival of the sound 
wave. The timing information from each receiver 210, 220 
and 230 is transmitted back to a central processing location 
250, where the timing information is processed to obtain the 
location of the lightning strike, see FIG. 3. The accuracy of 
the system can be enhanced in one embodiment by using 
more than three receivers, such as receiver 240. A network 
of four or more receivers, for example, can be used to 
resolve uncertainties introduced by wind speed. Echoes and 
reflections from objects within the monitored area can be 
removed using common digital signal processing tech- 
niques. 

The present invention allows for the accurate location of 
a lightning strike within a few’ meters using a combination 
of electric field and sonic sensors. One advantage of the 
present system includes the fact that high-speed digitizers 
are not required. The electric field pulse is used to start a 
time counter, and the sound wave is used to stop the time 
counter. A one-millisecond error in determination of the 
timing can result in an error of about 30 cm. Existing 
wideband, large baseline lightning location systems require 
timing accuracies better than a fraction of a nanosecond to 
achieve this kind of accuracy. The present invention pro- 
vides an inexpensive and easy to install system, with mini- 
mal maintenance and calibration requirements. 

The present invention allows for the determination of the 
distance to a lightning strike without requiring fast recording 
or digitizing equipment. Further, an algorithm to combine 
the information from a network of receivers allows for 
fine-tuning the system. Also, for example, an algorithm can 
be implemented to determine the location of a lightning 
strike when the wind is nonzero. This is important since 
nonzero wind will result in noncircular distance patterns 
around each receiver. 

In another embodiment, the present system discriminates 
between nearby and distant thunder. Because sound waves 
are attenuated as they propagate through air, with high 
frequencies decaying faster than low frequencies, the fre- 
quency spectrum of nearby thunder contains higher fre- 
quency components than thunder from a . distant lightning 
strike. The characteristic “rumble” from thunder consists 
mainly of frequencies below 100 Hz, while the “clap” 
occurred following a close lightning strike contains compo- 
nents above several kHz. 

The following receiver algorithm illustrates discrimina- 
tion between local and distant lightning strikes to avoid 
erroneous detection outside an area of interest: 


x =rf(“data) Read data file. 

h ^(“hpf 1 ) Read highpass filler, cutoff of 200 Hz, finite 

impulse response (FIR), 31st order. 

y =x$h Perform time-domain convolution between data file 

and filter impulse response. This is equivalent 
to multiplying the spectrum of the thunder data 
by the frequency response of the filter. The 
result of the convolution is stored in variable y. 

y =y-avg(y) Center variable y around zero (remove DC 

component). 

z =abs(y) Get the absolute value of y and storage in 

variable z. This is equivalent to performing a 
rectification of the filtered thunder signal. 

I =rf(“Ipf’) Read low pass filler, cut off of 10 Hz, finite 

impulse response (FIR), 31st order. 

z =z$l Perform a time-domain convolution between the 

rectified thunder signal and the low pass filter 
impulse response. This results in the envelope 
of the thunder signal. 
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-continued 


pop =FAI .SK 

Determine the start of U e thunder waveform by 
eliminating the 

5 

for(i=0:len: 1 > 

“'pop” sound caused by icarhy lightning. This 
is done by 

fori j:=0:500) 

eliminating pulses with \ duration shorter than 



50 ms. 


if(z|i+j|<thr) 
then pop=TRUT' 

end 


to 

end 

if(pop=FAt.SH) 

th under start=i 

The start of the thunder waveform is the amplitude 
*1" of the 


end 

data set. This is the vah e that is used to determine 



the 

15 

end 

spatial location of the lightning strike. 



The primary sources of location error in the described 
invention are due to variations in the speed of sound and due 
to the effect of wind. The sound speed in air is a function of 20 
the temperature (changing as the square root of the absolute 
temperature) and the molecular weight. The latter changes 
are due to variations in humidity, which can typically be 
ignored. By augmenting the present invention with exter- 
nally measured temperature, wind direction and wind speed, 25 
the accuracy of the system can be preserved. An environ- 
mental factors component 252 car be included with the 
present invention to input environmental data such as 
temperature, wind direction and wind speed. One embodi- 
ment of the described invention includes a temperature 
measurement at the central processor; which allows for the 
calculation of acoustic speed. In another embodiment, data 
is acquired from another source, such as a local meteoro- 
logical station to provide the input to compute sound speed. 
This computation is performed using widely known and M 
accepted equations. 

To correct for w ind speed, note that sound will travel with 
respect to the terrain as a sum of the wind speed and 
direction vector and the sound speed. Therefore the com- 
ponent of wind along the direction between the receiver and 4 
the lightning source can either retard or advance the effective 
velocity of the thunder. In one embodiment of the invention, 
three remote receivers are used and the central processor 
acquires external information on the wind speed and direc- 
tion from a local measuring station. The three receivers are 
used to compute a preliminary source location. The angles 
between the receivers and the estimated source location are 
used to recompute the estimated ranges, R, using the fol- 
lowing formula: 

R=TOA*V=TOA*ic-\\, cos (0-©)) 

where TOA=measured time interv il between the electro- 
magnetic and sonic signals 

V=effective sonic velocity 55 

c=actual sonic velocity based on temperature 
V M ,=wind speed 

0=azimuth direction measured from receiver to source 
©=azimuth of wind vector as measure from North 
These new range estimates are used to determine a new 60 
location as previously described. This process can be iter- 
ated until the estimated position of the source does not vary 
more than the expected variance based on the GDOR 
Alternatively, if wind speed and direction information are 
not available externally, the wind speed and direction can be 65 
estimated if a fourth (or more) receiver 240 is included. 
Again, the measured ranges are used to compute a source 


6 

location using iterative least squares procedure (i.e. guess a 
location, then use linearized range equations to derive a 
correction vector, etc.). Once this procedure has resulted in 
an estimate, a second, nonlinear programming technique 
operating on the range equation given above is used to 
estimate the location (x, y) of the source along with the wind 
velocity and direction. This procedure minimizes the vari- 
ance of all four variables jointly. This accuracy of this 
process improves as the number of receivers is increased. 
The efficacy of this nonlinear process is also improved by 
using the residuals (the differences between the range to the 
estimated point and the measured ranges) to estimate the 
wind speed and direction. 

The process is now described in detail. The equation for 
the range in terms of Xp and Yp the coordinates of the source 
point and Xi and Yi the coordinates of receiver “i” is given 
by (ignoring for the moment the wind): 

R. i 

it. = — = -<(A> - x,)- + o> - hr) 

c c 


wffiere there is one such equation for each detector. If this 
equation is expanded in terms of the coordinates of the event 
the following is obtained: 


( OR \ ( OR \ 

Ar ”= Arp+ (^P >+ l^P'' 


Take ATm to be the measured time difference while ATp 
is the time difference that would be measured from an 
assumed position p and evaluate the derivatives at the 
assumed position p, and interpret that the AX r and AY,> are 
components of a first order correction to the assumed 
position in a direction to reduce the difference between the 
measured ATm and the computed ATp. By listing these 
linear equations in rows a matrix equation can be con- 
structed and solve for the AX P and AY P values. The deriva- 
tive terms are: 


I OR \ _ Xp 
Ua>; r 

— 

\0 Yp) R 

which are easily calculated using the value for an assumed 
point Xp and Yp. Our total equation in matrix form looks 
like: 


' AT,; — ATp 1 

a To -Ar P = 
k AT.< - XTp , 


Xp 


Ap 

l 


Yp 


A’p Xp 

*7 *7 


AAp 

A)> , 


>p 

Rs 


In matrix notation this same equations reads: 
dXT-HAX 

To solve for the AX matrix the generalized inverse of the 
matrix H is taken. The solution is given by: 

AXMH T H)- ] H T bM 


Once the correction vector AX is solved, it can be added to 
the presumed values of Xp and Yp to create a new estimate. 
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Thus the estimate at time step n is transformed into a new 
estimate a step n+1. After each step, the size of the residuals, 
the elements of the 5AT vector, are evaluated. When they 
become sufficiently small the process can be considered 
complete. The final values of the Xp and Yp coordinates are 5 
used as the initial values of the next step which will estimate 
the event position jointly with wind speed and direction. 

Since each measurement of AT contains errors due to the 
wind speed, direction, and other errors (refraction, timing 
errors, errors in estimating the peak of the sound waveform, 10 
and so on) it is desirable to develop a figure of merit for the 
resulting solution. Specifically relating the errors in the AT's 
to the errors in and that represent our final and best 
estimates is desirable. The covariance matrix of the errors in 
the four AT measurements given as follows: 15 


This ellipse has semimajor and semiminor axes given by: 


a 2 = -(tr“ + tr;) + / -(ir? - it - x ) 2 +cr ^ 


b 1 = ^<<*7 +<rj ) - J - <r? ) 2 + <r^ v 

So unless the covariances are zero, the ellipse is inclined 
with respect to the x axis so that the angle between the 
semimajor axis and the x axis is given by: 

2<r u 

tan2y = — — — 

(r- - <r- 


comr.i = £{fiAT,4ni = 


I 


<r i 

<r, 2 


tr 14 

tr\2 

cr 5 


w 2 4 

<Hi 

{T$2 


4 

cr 4 , 

ti42 


<ri 


Note that the value of confidence applied was 68.3%, the 
probability that the actual value falls within the ellipse is 
significantly reduced to only 39.4%. Thus many practitio- 
ners derive an ellipse based on 2.447 times the one sigma 
values to get an ellipse encircling a 95% confidence level. 


where the diagonal terms represent the variances of each 
measurement and the off-diagonal terms represent the cova- 25 
riances between measurements. The off diagonal terms 
represent the degree of correlation between the measure- 
ments. In the case of random errors only (such a digitizing 
timing jitter and random variations in the sound waveform) 
these terms are zero, in the case of systematic errors 30 
including the effects of wind, they will not be zero. If the 
error sources are unknown, such as the wind speed and 
direction, it can be assumed that the covariances are zero. 
The matrix can be constructed by inserting identical values 
for the diagonal elements representing AT measurement 35 
errors derived from field tests or other considerations. The 
covariance matrix becomes the unit matrix multiplied by a 
scalar quantity oAT which represents the level of confidence 
chosen. Typically the sigma value is selected to represent a 
level of confidence such that the probability of a randomly 40 
chosen value of delta T falling within that distance of the 
mean is 68.3%. The following steps will result in error 
figures for the coordinates that will be based on the same 
confidence limit. The covariance matrix then is a unit matrix 
of rank four multiplied by a scalar value of error in units of 45 
time. 

Likewise, there is a covariance matrix of the estimated 
position values X r and Y^> given as follows: 


The GDOP is given by: 

GDOP = \lTRACE((H r H)-'H T ) 

This represents the magnification factor of the error based on 
the geometry of the lightning detectors and the computed 
lightning location. The GDOP is relevant to the present 
invention only in that should it be possible to select the best 
subset of available measurements (more than 4 receivers) 
the processor will select those that gave the best GDOP. 

A preliminary solution has been provided without mod- 
eling any wind effects. Next, estimates of the wind speed and 
direction from at least four lightning delta T measurements 
can be derived. By returning to the basic equation of delta 
T as a function of wind speed and direction: 


(c - VhsCOS(0 - <?)) 

Again, if a lightning location is assumed and a set of values 
for the wind speed and direction, this function can be 
expanded in a first order Taylor series as follows: 


COV(X) = E\e x e t x \ = ^ 

Y 


M ,r "j 

50 

1 dF v 

(OF \ 

( dF 

l <Tyx *T- } ) 


Ar ''= A7>+ (^} 

rfx+ U>) 

■Mssr. 


/V\v 


(£)-- 


where again, the diagonal terms are the variances of each 
coordinate and the off diagonal terms are their covariants ^ 
already derived transformation from the delta T equation to 
the delta X’s can be used to find the covariance matrix in X 
as follows: 

L x =(H T H)~ [ H Tl L ST 60 

Since the covariance matrix is symmetrical and positive 
definite it represents a quadratic form. If expressed in the 
following form it gives rise to an ellipse centered on the final 
estimates of X P and Y P : 65 

/(.V)=.V%Y 


where each derivative is evaluated at the assumed values of 
position and wind characteristics. If four such measurement 
equations are written in matrix form the following results: 


&T\ 


' 0 \ 

01 

y\ 

' 

/ dXp 

6 T 2 


ti 2 

02 

72 

f }2 

dY P 

st 


<*3 

0 i 

yy 


dV w 

,st 4j 


< ti4 

Pa 

y.i 

K - 

, d 0 , 


or ST = A AX 


where the A matrix elements are the values of the derivatives 
evaluated at the assumed values of the position, wind speed 
and direction. The equations for the derivatives are: 
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() 

R 

(X; - X P ) 

0 Xp 

|c(l - -^-cosftf - <h) J 

,R 

1 — COS(tf- (h) 


o: - r P ) 




r/?| 1 — cos(tf-rt)J 

RcosW-ii) 


£'"| 1 - -^-COS(W-rf)) 

«•-(* — ^-cos<0-tf)) 


3. The system of claim 1 wherein the processor compen- 
sates for wind speed and wind direction while determining 
the location of lighting strikes. 

4. The system of claim 3, wherein the processor compen- 
5 sates for wind speed and wind direction in the absence of 

wind speed and wind direction measurements. 

5. The system of claim 1 wherein the processor compen- 
sates for temperature while determining the location of 
lighting strikes. 

jo 6 . The system of claim 1 wherein the processor compen- 
* sates for humidity while determining the location of lighting 
strikes. 

7. The system of claim 1 wherein the plurality of receivers 
comprise three receivers. 

is 8 . A system to determine a location of lighting strikes 


When the values of the coefficients are evaluated, a simple 
matrix inversion is computed to ind the values of the 
elements of the AX vector. These four values are then added 
to the original set of assumed values for the four sought 
parameters to derive a new set. The residuals (elements of 
ST matrix) are evaluated between each iteration. When the 
residuals become sufficiently small, he process is complete. 
To compute the variances proceed us before with: 

V — A-Iv 

and the error ellipse and GDOP are computed as before with 
the exception of only using the upper left block of four 
elements of the sigma x matrix. 

This process results in optima) It ast square estimates of 
the position of the lightning strike and the errors in those 
estimates. 

Conclusion 


comprising: 

at least three receivers wherein each receiver comprises 
an electric field sensor, an acoustic sensor, and a 
controller to provide a receiver output indicating a 
20 calculated time differential between an electric field 
pulse and a sound wave; and 

a processor coupled to the receivers to determine the 
location of lighting strikes in response to the output 
from the receivers. 

- 9. The system of claim 8 wherein the processor compen- 

sates for environmental conditions including wind speed and 
wind direction, temperature and humidity while determining 
the location of lighting strikes. 

10. The system of claim 9, wherein wind speed is cor- 
™ rected for by 

using angles between the receivers and the lighting loca- 
tion to compute estimated ranges R, using the follow- 


A system and method of determining locations of light- 
ning strikes has been described. The system includes mul- ^ 
tiple receivers located around an area ot interest, such as a 
space center or airport. Each recei\er monitors both sound 
and electric fields. The detection of an electric field pulse 
and a sound wave are used to calculate a range circle around 
each receiver in which the lighting is detected. A processor 40 
is coupled to the receivers to accurately determine the 
location of the lighting strike. The processor can manipulate 
the receiver data to compensate for environmental variables 
such as wind, temperature, and humidity. Further, the system 
can discriminate between distant and local lightning strikes. 

Although specific embodiments have been illustrated and 
described herein, it will be appreciated by those of ordinary 
skill in the art that any arrangement, which is calculated to 
achieve the same purpose, may be substituted for the specific 
embodiment shown. This application is intended to cover 50 
any adaptations or variations of the present invention. 
Therefore, it is manifestly intended that this invention be 
limited only by the claims and the equivalents thereof. 

What is claimed is: 

1. A system to determine a location of lighting strikes 55 
comprising: 

a processor; and 

a plurality of receivers coupled t > the processor, wherein 
each of the plurality of receivers comprises an electric 
field sensor, an acoustic senior, and a controller to 60 
provide a receiver output indicating a calculated time 
differential between an electric field pulse and a sound 
wave, the processor determine s the location of lighting 
strikes in response to the output from the plurality of 
receivers. 

2. The system of claim 1 w'herei \ each of the plurality of 
receivers are located up to at least one kilometer apart. 


ing formula: 

R-TOA * V=TOA *( c — V w cos ( &-© ) 1 

where TOA=measured time interval between the electro- 
magnetic and sonic signals 
V=effective sonic velocity 
c=actual sonic velocity based on temperature 
V )r =wind speed 

0 =azimuth direction measured from receiver to source 
G=azimuth of wind vector. 

11. The system of claim 8 wherein each of the receivers 
are located up to one kilometer apart. 

12. A method for determining a location of lightning 
strikes comprising: 

locating a network of at least three electric field sensors 
and at least three sonic sensors in an area of interest; 
collecting lightning strike information, including a time of 
arrival of an electric field pulse and an associated sound 
wave from a lightning strike; 
processing the lightning strike information recorded by 
the at least three electric field sensors and the at least 
three sonic sensors, including measuring a time differ- 
ence between the arrival of the electric field pulse and 
the sound wave at each electric field and sonic sensor; 
and 

determining the location of the lightning strike. 

13. The method of claim 12 wherein determining the 
location of the lightning strike comprises compensating for 
environmental conditions including wind speed and wind 

65 direction, temperature and humidity. 

14. The method of claim 12 wherein the at least three 
electric field sensors are located up to one kilometer apart. 
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and the least three sonic sensors are located up to one 
kilometer apart. 

15. The method of claim 12 wherein determining the 
location of the lightning strike comprises comparing the 
arrival of the electric field pulse and the sound wave for each 5 
pair of receivers. 

16. The method of claim 12 wherein the network com- 
prises at least four electric field sensors and at least four 
associated sonic sensors. 

17. The method of claim 16, wherein wind speed is 10 
corrected for by 

using angles between the receivers and the lighting loca- 
tion to compute estimated ranges R, using the follow- 
ing formula: 

15 

R=TOA*V=TOA*U-V\, cos (6-0)) 

where TOA=measured time interval between the electro- 
magnetic and sonic signals 
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V=eflFective sonic velocity 
c=actua! sonic velocity based on temperature 
V M ,=wind speed 

0=azimuth direction measured from receiver to source 
0=azimuth of wind vector. 

18. The method of claim 12 wherein the processor com- 
pensates for wind speed and wind direction while determin- 
ing the location of lighting strikes. 

19. The method of claim 12, wherein the processor 
compensates for wind speed and wind direction in the 
absence of wind speed and wind direction measurements. 

20. The method of claim 12 wherein processing the 
lightning strike information comprises discriminating 
between lightning strikes based upon distance to filter out 
lightning strikes outside the area of interest. 

* * * # * 
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